NASA  Contractor  194901 
ICASE  Report  No.  94-23 


AD-A280  580 

HBHIRlilff 


m  ICASE 


A  PHASE  EQUATION  APPROACH  TO 
BOUNDARY  LAYER  INSTABILITY  THEORY: 
TOLLMIEN  SCHLICHTING  WAVES 


ptf*  document  he* _ 

tw  public  and  ialac  in 


Contract  NAS1-19480 
April  1994 


mc  quality  inspected^ 

94-19273 

lllllllll 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  VA  23681-0001 


Operated  by  Universities  Space  Research  Association 

if*  6  23 


089 


ICASE  Fluid  Mechanics 


Due  to  increasing  research  being  conducted  at  ICASE  in  the  field  of  fluid  mechanics, 
future  ICASE  reports  in  this  area  of  research  will  be  printed  with  a  green  cover.  Applied 
and  numerical  mathematics  reports  will  have  the  familiar  blue  cover,  while  computer  science 
reports  will  have  yellow  covers.  In  all  other  aspects  the  reports  will  remain  the  same;  in 
particular,  they  will  continue  to  be  submitted  to  the  appropriate  journals  or  conferences  for 
formal  publication. 


Accesfon  For  \ 

. . - . 

- J 

NTIS  CRA&l  V/ 

Dll  C  7  A3 

! 

Ur.ji;::ou.:ced 

j 

Justification 

f 

By 

• 

I 

Distribution  / 

. 1 

* 

Availability  Cooes  *1 

A  PHASE  EQUATION  APPROACH  TO 
BOUNDARY  LAYER  INSTABILITY  THEORY: 
TOLLMIEN  SCHLICHTING  WAVES 


Philip  Hall1 

Department  of  Mathematics 
University  of  Manchester 
Manchester  M13  9PL,  U.K 


Abstract 

Our  concern  is  with  the  evolution  of  large  amplitude  Tollmien-Schlichting  waves  in 
boundary  layer  flows.  In  fact  the  disturbances  we  consider  are  of  a  comparable  size  to 
the  unperturbed  state.  We  shall  describe  two-dimensional  disturbances  which  are  locally 
periodic  in  time  and  space.  This  is  achieved  using  a  phase  equation  approach  of  the  type 
discussed  by  Howard  and  Kopell  (1977)  in  the  context  of  reaction-diffusion  equations.  We 
shall  consider  both  large  and  0(1)  Reynolds  numbers  flows  though,  in  order  to  keep  our 
asymptotics  respectable,  our  finite  Reynolds  number  calculation  will  be  carried  out  for  the 
asymptotic  suction  flow.  Our  large  Reynolds  number  analysis,  though  carried  out  for  Bla- 
sius  flow,  is  valid  for  any  steady  two-dimensional  boundary  layer.  In  both  cases  the  phase 
equation  approach  shows  that  the  wavenumber  and  frequency  will  develop  shocks  or  other 
discontinuities  as  the  disturbance  evolves.  As  a  special  case  we  consider  the  evolution  of  con¬ 
stant  frequency/ wavenumber  disturbances  and  show  that  their  modulational  instability  is 
controlled  by  Burgers  equation  at  finite  Reynolds  number  and  l>y  a  new  integro-differential 
evolution  equation  at  large  Reynolds  numbers.  For  the  large  Reynolds  number  case  the 
evolution  equation  points  to  the  development  of  a  spatially  localized  singularity  at  a  finite 
time.  The  three-dimensional  generalizations  of  the  evolution  equations  is  also  given  for  the 
case  of  weak  span  wise  modulations. 
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1  Introduction 


Most  boundary  layers  of  practical  importance  are  susceptible  to  a  variety  of  instability  mech¬ 
anisms  which  lead  to  the  onset  of  transition  to  turbulence.  Usually  more  than  one  mechanism 
will  be  operational  in  any  particular  case,  and  a  full  understanding  of  how  transition  occurs  will 
require  a  detailed  understanding  of  the  nonlinear  interaction  of  the  different  modes  of  insta¬ 
bility.  Here  we  will  concern  ourselves  with  the  strongly  nonlinear  evolution  of  a  slowly- varying 
Tolbnien-Schlichting  wave  system.  In  the  first  instance  we  consider  lower  branch  Tollmien- 
Schlichting  waves  which  are  known  to  be  governed  by  triple-deck  theory  e.g.  Smith  ( 1979a, b), 
Hall  and  Smith  (1984),  Smith  and  Burggraf  (1985).  Then  we  shall  consider  the  corresponding 
problem  at  finite  Reynolds  numbers.  It  is  worth  pointing  out  that  the  approach  we  use  here 
can  be  used  to  describe  other  modes  of  instability  and  is  in  fact  based  on  ideas  given  some 
years  ago  by  Howard  and  Kopell  (1977)  who  were  interested  in  the  evolution  of  nonlinear  wave 
systems  in  reaction- diffusion  equations. 

The  first  application  of  triple-deck  theory  to  describe  the  linear  and  nonlinear  growth 
of  lower  branch  Tollmien-Schlichting  waves  is  apparently  due  to  Smith  (1979a, b)  though 
Lin  (1966)  clearly  recognized  the  appropriate  large  Reynolds  number  scalings  for  Tollmien- 
Schlichting  waves  long  before  triple-deck  theory  was  invented.  The  investigation  of  Smith 
(1979a)  showed  how  nonparallel  effects  could  be  taken  care  of  in  a  self-consistent  manner  us¬ 
ing  asymptotic  methods.  Previously  Gaster  (1974)  used  a  successive  approximation  procedure 
to  tackle  the  same  kind  of  problem.  Subsequently  Smith  (1979b)  showed  how  the  nonlinear 
growth  of  Tollmien-Schlichting  waves  could  be  taken  care  of  using  triple-deck  theory.  How¬ 
ever  the  results  of  Smith  (1979b),  and  the  subsequent  extension  to  three-dimensional  modes 
by  Hail  and  Smith  (1984),  are  confined  to  the  weakly  nonlinear  stage  where  an  ordinary  dif¬ 
ferential  amplitude  equation  describes  the  initial  stage  of  the  bifurcation  from  a  disturbance 
free  state.  Some  years  later  Smith  and  Burggraf  (1985)  discussed  the  high  frequency  limit  of 
the  lower  branch  triple-deck  problem  and  uncovered  a  sequence  of  nonlinear  structures  gov¬ 
erning  a  sequence  of  successively  more  nonlinear  wave  interactions.  Subsequently  Smith  and 
Stewart  (1987)  investigated  the  interaction  of  three-dimensional  modes  at  high  frequencies 
and  obtained  excellent  agreement  with  the  experiments  of  Kachanov  and  Levechenko  (1984), 
though  Khokhlov  (1994)  in  a  recent  theoretical  investigation  argues  that  the  work  of  Smith 
and  Stewart  (1987)  needs  some  modification. 

In  the  first  instance  we  shall  restrict  our  attention  to  two-dimensional  waves  and  determine 
how  the  wavenumber  and  frequency  of  a  wavesystem  may  be  found  as  it  moves  through  a 
growing  boundary  layer.  This  problem  has  not  yet  been  addressed.  Intuitively  one  would 
expect  that  a  small  amplitude  wave  would  evolve  from  its  weakly  nonlinear  form  into  a  larger 
amplitude  state  until  it  is  described  by  the  Smith-Burggraf  structure  at  sufficiently  large  values 
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of  the  local  frequency  of  the  disturbance.  Our  calculations  show  that  this  is  not  the  case  and 
indeed  show  that  at  high  frequencies  locally  periodic  forms  of  the  modes  found  by  Smith  and 
Burggraf  probably  do  not  exist.  Certainly  it  would  appear  that  they  do  not  connect  to  the 
weakly  nonlinear  state  discussed  by  Smith  (1979b). 

The  asymptotic  structure  we  use  is  based  on  the  so-called  ‘phase-equation  ’approach  used 

so  successfully  to  describe  large  amplitude  Benard  convection  in  large  containers  by,  amongst 

others,  Kramer,  Ben  Jacob,  Brand  and  Cross  (1982),  Cross  and  Newell  (1984),  Newell,  Passot 

and  Lega  (1993).  Using  this  approach  it  has  been  possible  to  describe  the  experimentally 

observed  slowly  varying  planform  of  Benard  convection.  Thus,  for  example,  the  dislocation 

of  convection  rolls  is  now  reasonably  well  understood  using  the  phase  equation  approach. 

% 

Interestingly  enough  it  turns  out  that  the  essential  ideas  of  this  approach  had  been  elucidated 
in  the  context  of  traveling  wave  instabilities  several  years  earlier,  see  Howard  and  Kopell  (1977) 
and  indeed  the  method  can  be  found  in  Whitham  (1974).  The  evolution  of  traveling  waves  in 
a  Blasius  boundary  layers  is  the  subject  of  the  first  part  of  this  work  and  not  surprisingly  the 
analysis  to  be  used  has  similarities  with  that  of  the  latter  authors. 

The  essential  idea  behind  the  phase  equation  approach  may  be  explained  in  the  following 
manner.  Suppose  there  exists  some  flow  which  is  unstable  to  a  traveling  wave  disturbance  of 
wavenumber  a  and  frequency  ft.  For  a  fully  nonlinear  disturbance  the  frequency  ft  will  be 
a  function  of  a  which  itself  can  be  thought  of  as  a  function  of  A,  a  measure  of  the  size  of 
disturbance.  If  we  let  A  tend  to  zero  then,  for  small  A,  the  quantities  a  and  ft  will  differ  from 
their  linear  neutral  values  by  0(A)2  so  that  finite  amplitude  disturbances  begin  as  supercritical 
bifurcations  from  the  basic  state.  For  0(1)  values  of  A  the  quantities  a  and  ft  are  accessible 
only  by  numerical  means,  see  for  example  Herbert  (1977)  for  details  of  the  computation  of  a 
and  ft  for  Tollmien-Schlichting  waves  in  plane  Poisseuille  flow  or  Conlisk,  Burggraf  and  Smith 
(1987)  for  a  similar  calculation  for  Tollmien-Schlichting  waves  in  Blasius  boundary  layers  at 
large  Reynolds  numbers.  In  some  cases  the  frequency  of  the  waves  is  zero  and,  the  wavenumber 
of  the  disturbances  may  be  sensibly  held  fixed  when  the  control  parameter  or  disturbance  size 
is  varied,  see  for  example  Hall  (1988)  for  a  discussion  of  the  fully  nonlinear  Gortler  problem  in 
a  growing  boundary  layer.  For  a  traveling  wave  disturbance  in  a  growing  boundary  layer  we 
expect  that  the  wavenumber  and  frequency  of  the  disturbance  should  change  as  it  propagates 
into  locally  less  or  more  unstable  parts  of  the  flow.  The  phase  equation  approach  provides 
a  rational  framework  for  following  such  an  evolution.  If  the  wave  has  local  frequency  and 
wavenumber  which  are  0(1)  with  respect  to  the  variables  t  and  x  we  introduce  slowly  varying 
variables  T  and  X  by  writing 

T  =  St,  X  =  6x, 

and  we  now  think  of  a  and  ft  as  being  functions  of  X  and  T.  Thus  we  may  introduce  the 
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phase  function  Q(X,T)  defined  by  0  =  $(X,T)/S  with 

a  =  &x ,  ft  =  —Ox 

and  as  a  consequence  the  wave  system  evolves  such  that 


c*t  +  =  0. 


(1.1) 


Partial  derivatives  with  respect  to  x  and  t  must  then  be  replaced  using 

d_ 
dx 


3  d  d  d  d 

ade+6dx'  di  ~nd®  +  6df' 


and  if  we  then  equate  terms  of  0(6)°  we  recover  the  unmodulated  equations  of  motion  with  a 
and  ft  playing  the  role  of  wavenumber  and  frequency.  Thus  the  leading  order  problem  using 
the  phase  equation  approach  is  simply  the  unmodulated  case  with  a  and  ft  being  functionally 
related  in  order  that  the  system,  with  a  given  disturbance  size,  has  a  solution.  At  next  order 
(usually  0(6)  but  in  fact  0(61/3)  in  triple  deck  problems)  we  obtain  a  linearized  inhomogeneous 
form  of  the  leading  order  problem.  Due  to  the  invariance  of  the  problem  under  a  translation 
in  the  x  direction  it  is  easy  to  see  that  the  linearized  homogeneous  form  of  this  system  has  a 
non-trivial  solution  so  that  some  solvability  condition  must  be  satisfied  if  the  inhomogeneous 
problem  is  to  have  a  solution.  This  solvability  condition  is  satisfied  by  introducing  an  expansion 
of  ft  in  appropriate  powers  of  6.  The  solution  of  this  problem  then  enables  us  to  write  down  the 
asymptotic  form  of  (1.1)  upto  the  second  order.  This  procedure  can  be  continued  in  principle 
to  any  order  and  the  coefficients  in  the  expansion  of  ft  are  found  as  solvability  conditions  at 
each  order.  The  evolution  of  a  given  wave  system  can  then  be  found  by  the  solution  of  the 
calculated  asymptotic  approximation  to  (1.1).  We  find  that  (1.1)  takes  on  a  particularly  simple 
form  if  the  wave  system  has  fixed  wavenumber  and  frequency  at  leading  order.  We  shall  in  this 
paper  calculate  (1.1)  correct  upto  second  order  for  both  wavesystems  governed  by  triple-deck 
theory  and  those  satisfying  the  two-dimensional  Navier-Stokes  equations.  We  can  show  that 
at  high  Reynolds  numbers  (1.1)  may  then  be  reduced  to  the  form 


dA  dA  d  [ c 
or  d£ 


d\ 

ds 


rds. 


(1-2) 


dtJ<  (*-01/3 

This  is  in  effect  the  evolution  equation  for  a  wavepacket  of  large  amplitude  Tollmien 
Schlichting  waves  in  a  growing  boundary  layer. 

At  finite  Reynolds  numbers  the  modulation  equation  corresponding  to  (1.2)  is 


At  +  AA^  =  iA^. 


(1-3) 


Thus  at  finite  Reynolds  numbers  Burgers  equations  controls  the  slow  dynamics  of  a  two- 
dimensional  wavesystem.  We  will  show  that  (1.2)  and  (1.3)  can  be  generalized  to  allow  for  a 
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weak  spanwise  dependence  of  the  modulation  and,  surprisingly,  it  turns  out  that  the  general¬ 
izations  in  both  cases  are  achieved  by  adding  a  term  proportional  to  to  both  equations. 
Here  £  is  a  slow  spanwise  variable. 

The  procedure  adopted  in  the  rest  of  this  paper  is  as  follows:  In  §2  we  derive  the  phase 
equation  for  two-dimensional  triple-deck  problems.  In  §3  we  describe  the  numerical  work 

i 

required  to  determine  the  quantities  appearing  in  the  equation.  In  §4  we  look  at  the  special 
case  of  almost  uniform  wavetrains  and  derive  (1.2).  In  §5  we  show  how  (1.2)  can  be  derived 
by  a  more  conventional  multiple  scale  approach  directly  from  the  triple-deck  equations.  The 
phase  equation  approach  for  a  boundary  layer  at  finite  Reynolds  numbers  is  then  discussed 
in  §6.  The  modulation  equation  (1.3)  is  derived  in  that  section  as  a  special  case  for  almost 
uniform  wavetrains.  Finally  in  §7  we  draw  some  conclusions  and  give  the  generalized  form  of 
the  evolution  equations  which  account  for  weak  spanwise  modulations  and  nonparallel  effects. 

2  Derivation  of  the  phase  equation  for  2D  triple  deck  prob¬ 
lems 

Our  concern  is  with  the  structure  of  fully  nonlinear  solutions  of  the  triple-deck  equations  gov¬ 
erning  the  evolution  of  two-dimensional  Tollmien-Schlichting  waves  in  incompressible  boundary 
layers.  Following  the  usual  notation,  e.g.  Smith  (1979a),  the  appropriate  differential  equations 
in  scaled  form  are 


du  dv 
dx  dy  ’ 

(2.1a) 

du  du  du  dp  d2u 

dt  U  dx  +  V  dy  dx  +  dy2  ’ 

(2.1b) 

which  must  be  solved  subject  to 

u  =  v  =  0,  y  =  0, 

(2.2a) 

u  ~  y  +  A(x,t),  y  -»•  oo, 

(2.2b) 

1  roo  &A 

P=  4  (  ds  ,ds. 

i r  j-o o  (a:  -  s) 

(2.2c) 

The  displacement  function  A  and  the  pressure  p  depend  only  on  x  and  t  and  if  we  wish  to 
consider  other  flows  the  pressure  displacement  law  (2.2c.)  must  be  modified  accordingly.  Linear 
Tollmien-Schlichting  waves  correspond  to  perturbing  u  in  the  form 


u  =  y+U(y)eiaix-%},  (2.3) 


with  U  small  and,  from  Smith  (1979a),  the  eigenrelation  takes  the  form 

too 

K  (&)  =  (ia)1/3a  /  A^dr). 


(2.4) 
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where  A,  is  the  Airy  function  and  £0  =  -7 — - 

(ta 

real  show  that  neutral  stability  occurs  for  f l  ~  2.298,  a  ~  1.001,  and  that  ft*  is  positive  for  all 
frequencies  greater  than  the  neutral  value.  In  the  high  frequency  limit  it  can  be  shown  from 


)2/3’ 


Solutions  of  (2.4)  with  ft  complex  and  ft 


(2.4)  that 


ft  =  ft  +  0(ft~1/2) 


(2.5) 


The  above  limit  was  discussed  in  detail  by  Smith  and  Buggraff  (1985)  who  investigated  the 
possible  nonlinear  structures  which  emerge  in  that  limit.  The  structures  found  by  Smith  and 
Burggraf  depend  crucially  on  the  fact  that  the  right  hand  side  of  (2.5)  is  complex  only  at  order 
ft-1/2  so  that,  even  though  a  wave  is  never  neutral,  its  growth  can  be  balanced  at  higher  order 
by  nonlinear  effects.  Here  our  interest  is  with  the  case  a  =  0(1),  ft  =  0(1),  but  we  shall  allow 
for  a  slow  evolution  of  the  wavesystem  as  it  moves  through  the  boundary  layer.  The  essential 
details  of  our  approach  are  to  be  found  in  Howard  and  Kopell  (1977)  who  were  concerned  with 
slowly  varying  waves  in  reaction  diffusion  systems.  As  a  first  step  we  introduce  slow  time  and 
space  variables,  T  and  X,  by  writing 

T  =  St,  (2.6a) 

X  =  St,  (2.6b) 


where  6  is  a  small  positive  parameter.  We  shall  investigate  the  evolution  of  a  fully  nonlinear 
wavelike  solution  of  (2.1)-(2.2),  but  allow  the  wavenumber  and  frequency  to  be  slow  functions 
of  X  and  T. 

In  order  to  describe  such  a  structure  we  introduce  a  phase  function  0(x,  t)  such  that  the 
wavenumber  and  frequency  of  the  wave  are  defined  by 


d0  n  00 

a  ~  dx'  dt * 

The  wavenumber  and  frequency  must  therefore  satisfy 

da  dft_ 
dt  +  dx  ’ 


(2.7) 


(2.8) 


and  (2.8)  therefore  corresponds  to  the  conservation  of  phase.  Now  we  shall  assume,  following 
Howard  and  Kopell  (1977),  that  ft  and  ft  are  functions  only  of  X  and  T.  In  that  case  (2.8) 
reduces  to 

— +  — =  0  (29) 

dT  +  dx  1  j 


and  it  is  then  convenient  to  write  to  phase  variable  0  =  6  l0(X,T).  The  x  and  t  derivatives 
in  (2.1)-(2.2)  must  then  be  transformed  according  to 


d  3  .  d 

dx  ade  +6dx  ’ 

9  _  d  sd 

dt adQ+SdT ■ 


(2.10a) 

(2.10b) 
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We  seek  a  locally  wavelike  solution  of  (2.1)-(2.2)  and  impose  periodicity  in  the  phase  variable 
0.  It  remains  now  for  us  to  find  a  small  6  solution  of  the  full  triple  deck  problem  (2.1  )-(2.2). 
At  first  sight,  in  view  of  (2.10),  we  would  expect  to  develop  a  solution  of  that  system  in  terms 
of  6.  However,  it  turns  out  that  for  6  <<  1  the  leading  order  approximation  to  (2.1  )-(2.2)  has 
a  mean  term  correction  which  depends  on  the  slow  variable  X .  This  mean  flow  is  reduced  to 
zero  in  an  outer  0(£-1/3)  boundary  layer.  The  expansions  must  therefore  proceed  in  powers 


of  61/3  and  we  therefore  write 

ft  —  ft0  +  S  ^3ftj  -f-  •  •  • ,  (2.11a) 

u  =  u0  (X,T,Q,Y)  +  6''3ui(X,T,Q,Y)  +  ---,  (2.11b) 

v  =  uo(X,T,0,y)  +  61/3u1(X,T,0,T)  +  .--,  (2.11c) 

p  =  p0(X,T,Q)  +  61'3p1(X,T,  0)+-..,  (2.1  Id) 

A  =  Aq(X,T,  0)  +  61/3A\(X,T,  ©)+•••.  (2.1  le) 


The  leading  order  problem  is  then  found  to  be 


du0  dv0 

+  -5—  =  0, 
a0  ay 

(2.12a) 

_  du0  du0  du0  dpo  d'2u0 

n“ae+““"«e+^a#  =  °»e+v’ 

(2.12b) 

Uo  =  V0  =  0,  y  =  0, 

(2.12c) 

«o  ~  V  +  A0(X,  Q,T)  y  -*■  00, 

(2.1 2d) 

1 

Po  ax  /Loo  (0  -  y>y 

(2.12e) 

Hence  the  leading  order  problem  is  obtained  from  the  full  two-dimensional  problem  by  restrict¬ 
ing  attention  to  solutions  in  the  form  of  traveling  waves  of  local  wavenumber  a  and  frequency 
ft0.  This  specifies  a  nonlinear  eigenvalue  problem 


fto  —  flo  ( ^— )  ? 


(2.13) 


which  must  be  determined  numerically.  At  this  stage  we  assume  that  (2.13)  and  the  corre¬ 
sponding  nonlinear  eigenfunctions  uo,vo,po  and  Ao  are  known.  We  further  note  Aq  may  be 
written  in  the  form 

A0  =  Ao(X,T)  +  A0(Q,X,T)  (2.14) 


where  Ao  has  zero  mean  with  respect  to  0.  It  is  necessary  to  reduce  Ao  to  zero  before  the  main 
deck  is  encountered,  therefore  an  outer  boundary  layer  is  required.  Before  we  investigate  the 
outer  boundary  layer  it  is  convenient  to  discuss  the  next  order  system  in  the  y  =  0(1)  region. 


The  equations  to  be  satisfied  are 


dv\ 

dy 


=  0, 


(2.15a) 
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O  ^Ul  4. 

~svae  + 


,“°5e +“' 


subject  to 


dui  duo  dpi  d2ui  du0 

+  V°  dy  +Vl  dy  +  a  dO  dy 2  dQ  ’  ^2‘15  * 


U\  =  v\  =  0,  y  =  0,  (2.15c) 


pi  =  —  f  r^ds-  (2-isd) 

airj-oo  (x  -  s) 

In  addition  we  require  a  condition  involving  uj  at  the  edge  of  the  boundary  layer.  On  the  basis 
of  (2.2b)  we  might  expect  that,  u\  -*  Ai,  y  —*  oo,  is  the  appropriate  condition.  However  Ai  is 
essentially  the  displacement  function  in  the  main  deck  and  uj  is  modified  in  the  outer  boundary 
layer  in  which  the  mean  flow  correction  on  the  slow  scale  is  reduced  to  zero.  Therefore  we 
must  write  down  the  condition 


-*  Du  V  °°, 


(2.16) 


where  D\  is  to  be  found  in  terms  of  A\  by  a  consideration  of  the  outer  boundary  layer.  For 
large  values  of  y  we  have  u  ~  y  so  that  the  thickness  of  the  outer  layer  is  fixed  by  the  balance 


Hence  we  write 


dX  ~  dy2' 


r)  =  S1/3y 


and  now  develop  an  asymptotic  solution  of  (2.1)-(2.2)  valid  for  t?  =  0(1).  The  solution  here  is 
similar  to  that  in  the  main  deck  for  two-dimensional  triple-deck  problems.  We  write 

u  —  +  um(X,ti,T)  +  •  •  •}  +  jc^o  +  +  •  •  •} , 

«  =  {62/3vM(X,rhT)+--)  +  -^{Vo  +  61'3Vl+  ■■}. 

Here  the  terms  in  the  first  bracket  of  each  expansion  correspond  to  the  mean  flow.  If  we 
substitute  the  above  expansions  into  (2.1)-(2.2),  and  make  the  appropriate  changes  of  variables, 
then  we  find  that  the  function  Uo  is  given  by 

U0  —  Aq  +  6^!3H\  +  61/3AoUmv  +  •  •  • 


where  Hi  is  to  be  determined.  The  functions  um  and  dm  are  found  to  satisfy 


d2u  M 
dr/2 

9um 

dX 


duM 

-  nHx  M  =  °’ 


7 


which  are  to  be  solved  subject  to 


UM  =  Aq,  vm  =  0, 7/  =  0,  uM  *  0, 7/  *  oc. 

The  solution  of  the  above  problem  for  the  mean  flow  correction  is  most  easily  obtained  using 
a  Fourier  transform  with  respect  to  X .  We  find  that  the  function  um  is  such  that 


=  0)  =  J1/  ,-  / 
r  (s)  y“ 


X  §Ai 


oo  (s  -  X)z 


(2.17) 


If  we  now  match  the  0(61/3)  correction  to  the  wavelike  parts  of  the  expansion  between  the 
y  =  0(1),  y  =  0(^/3)  regions  we  obtain 


3iAo  fx  ^  .  D 

fl,  +  F(| )L 


whilst  letting  j)  -*  oo  in  the  definition  of  U\  gives 


tfi  =At. 


(2.18) 


This  closes  the  problem  for  the  order  6 1/3  problem  in  the  lower  layer.  Thus  in  addition  to 
(2.14),  (2.15)  and  (2.16)  we  have 


Di  =  Ax  + 


3  s  Aq 


p  i—xy/* 


(2.19) 


It  should  also  be  pointed  out  that  the  condition  that  the  mean  flow  correction  um  goes  to  zero 
at  infinity  can  be  relaxed,  if  required  for  some  Ao,  to  the  condition  that  um  tends  to  a  constant 
at  infinity.  This  does  not  alter  the  condition  (2.19).  The  quantity  ftj  is  now  determined  as  a 
solvability  condition  on  the  inhomogeneous  system  specified  by  (2.15),  (2.16)  and  (2.19).  Such 
a  condition  is  required  because  of  the  translational  invariance  of  the  leading  order  wavesystem 
with  respect  to  0.  Therefore  a  solution  of  the  homogeneous  problem  is  found  by  setting 
(«i,  V\,p\,  A\)  =  ^^(tiojUojPo^o).  In  order  to  find  the  appropriate  solvability  condition  it 
is  convenient  to  define  Z  =  (pi,ui,«i,«iy)T.  We  then  must  determine  the  condition  that  the 
inhomogeneous  system  given  below  has  a  solution: 

f  -  flZ+c,I  +  n>F<- 


B  = 


0  0 
0  0 


0  tioy  ««oe  v0 
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■  0 

0 

0 

0  ' 

0 

0 

0 

-a 

0 

,F,  = 

0 

■*“ 

0 

0 

0 

0 

0 

a 

0 

otio  —  fto 

0  . 

.  U0B  . 

subject  to 


u\  =  uj  =  0,  y  =  0, 


«1  -+  + 


3  rA _ * 

2  Joo  ( S  — 


*  ^ 


P  (§)  -Zoo  -  X)1/3 

1  foo 

Pi  =  -  /  li~ 

07T  J_ 


ds,  y  -*  oo, 


-00  *  ^1 

The  system  adjoint  to  the  homogeneous  form  of  the  above  problem  is 

-  DJ  -  CT  — 
dy  ~  UJ  1  00’ 


with 

‘  0  0  0  0 
_  0  0  0  «0y 

0  0  0  0 
0  0  1  uo 
J  =  (M,N,S,T)t, 


and  subject  to  the  conditions 

M  =  T  =  0,  y  =  0, 

S  —>  0,  y-*  oo, 

M  -  A/oo(0),  T  -  TooV*2,  y  —  oo  with  7^  =  (x~s)ds' 

The  condition  that  the  problem  for  {u\,v\,p\,A\)  has  a  solution  is  then  found  to  be 


with 

At  this  stage  we  can 
obtain 


fli 


rA  a" 

-kmLi2 


X  dAo  da 

- ds, 


3^ 


K(a)  = 


(s-xy/ 3 

Jo'  Moo(0)pod& 


(2.20) 


(2.21) 


Cfo°°T%£dedy 

write  down  the  phase  conservation  equation  correct  to  order  S'1/3.  We 


(2.22) 


and  the  expansion  procedure  given  above  can  in  principle  be  continued  to  any  order.  We 
postpone  a  discussion  of  the  implications  of  (2.22)  until  we  have  described  the  results  of  the 
calculations  required  to  determine  a, Ho, Hi. 
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3  The  numerical  work 


The  system  (2.12)  is  periodic  in  0  so  we  seek  a  solution  by  expanding  for  example  u0  in  the 
form 

OO 

vo  =  H  v0m(y)e,me  (3.1) 

—  OO 

After  eliminating  pgm  and  some  linear  terms  proportional  to  uo,n  from  the  0  momentum 
equation,  the  equation  to  determine  vom  may  be  written  in  the  form 

d4v0m 


dy4 


ilTl  (qiUqo  fl)  t'Om  —  Fyn  , 


(3.2) 


(3.3) 


where  Rrtl  is  a  nonlinear  function  of  {uom},  {ro,n}.  The  equation  for  the  mean  part  of  uq  is 
then  written  as 

d2ugg  _  dl 

dy 2  dy ’ 

where  /  is  a  nonlinear  function  of  {u0m},  { vgm } .  We  used  central  differences  to  evaluate  the 
derivatives  in  (3.2),  (3.3)  and  a  solution  of  the  resulting  nonlinear  system  was  found  by  iteration 
after  first  restricting  m  to  be  less  than  say  M.  In  our  iterative  technique  the  right  hand  sides 
of  (3.2)  and  (3.3)  were  evaluated  at  the  previous  level  of  iteration  and  one  boundary  condition 
was  replaced  by 

Uq(oo)  =  Co  +  ido 


where  c0  and  d0  are  prescribed  real  constants.  After  the  iterations  converged  we  then  adjusted 
a  and  fl  until  the  previously  ignored  boundary  condition  was  satisfied.  Note  here  that  ,  because 
the  solution  of  (2.12)  is  unique  only  upto  a  phase  shift,  the  values  of  a  and  fl  obtained  by  this 
procedure  are  functions  only  of  {cq  +  dg}1^2. 

The  grid  size  and  the  value  of  M  were  varied  until  converged  results  were  obtained.  In  the 
following  discussion  the  results  presented  correspond  to  M  =  32,  and  200  grid  points  in  the  y 
direction  with  “infinity”  at  y  =  10.  In  Figures  1  and  2  we  show  the  dependence  of  a  and  fto  on 
the  quantity  \Jcg  +  dg  which  is  a  measure  of  the  disturbance  size.  We  see  that,  as  predicted  by 
Smith  (1979b),  finite  amplitude  motion  begins  as  a  supercritical  bifurcation  from  Blasius  flow. 
In  Figure  3  we  plot  fig  as  a  function  of  a.  The  calculations  could  not  be  continued  beyond  the 
point  F  shown  on  the  Figure  ;  we  will  return  to  this  point  later.  We  further  notice  that  flo  is  a 
multiple  valued  function  of  a  for  a  range  of  values  of  a  and  that  fio(tv)  becomes  infinite  when 
a  ~  1.0145.  In  Figure  4  we  show  the  shear  stress  as  a  function  of  0  for  a  range  values  of  flo- 
The  results  shown  in  this  picture  suggest  a  reason  why  Figure  3  cannot  be  continued  beyond 
the  point  F.  We  see  that  as  the  point  F  is  approached  the  shear  stress  approaches  zero  at  a 
point.  In  Figure  5  we  show  how  the  contribution  to  the  shear  stress  from  the  different  modes 
varies  as  fio  varies.  We  see  that  the  the  higher  modes  grow  rapidly  as  F  is  approached.  This 
suggests  that  the  shear  stress  becomes  singular  as  F  is  approached. 
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Beyond  the  point  F  our  calculations  failed  to  converge  because  the  procedure  used  to  solve 
(3.2)-(3.3)  failed  to  drive  the  residuals  below  the  tolerance  level,  1 0 — 1 2 ,  used  throughout  our 
work.  A  similar  result  was  found  by  Conlisk  et  al  (1987)  who  solved  (2.12)  by  an  indirect 
method.  In  their  calculation  the  Tollmien-Sclichting  waves  were  first  forced  by  a  wall  motion 
and  then  their  properties  extrapolated  as  the  forcing  was  reduced  to  zero.  The  results  of 
Conlisk  et  al  (1987)  have  been  plotted  in  Figure  3  and  we  see  that,  on  the  whole,  there  is 
good  agreement  with  our  work  when  our  program  produced  converged  results.  Some  of  the 
larger  amplitude  results  by  Conlisk  et  al  were  obtained  by  reducing  the  tolerance  level  in 
their  iteration  procedure.  A  similar  reduction  of  the  tolerance  level  in  our  code  enables  us  to 
continue  Figure  3  for  slightly  larger  disturbance  amplitudes  but  we  do  not  plot  them  for  the 
following  reasons.  Firstly  we  found  that  a  reduction  of  the  tolerance  level  made  our  results 
very  sensitive  to  the  grid  size.  Secondly  a  reduction  of  the  tolerance  level  at  best  only  enabled 
us  to  continue  our  calculation  until  a  was  reduced  to  about  1.018.  Furthermore  the  results  of 
Figure  4  suggest  to  us  that  the  curve  of  Figure  3  terminates  at  a  point  close  to  F  where  all 
the  harmonics  are  excited  and  a  singularity  has  been  encountered.  Therefore  it  does  not  seem 
sensible  to  plot  results  obtained  by  reducing  the  tolerance  level  further.  Further  calculations 
were  carried  out  at  large  frequencies  in  order  to  find  finite  amplitude  solutions  of  the  type 
predicted  by  the  Smith- Burggraf  theory.  Despite  a  careful  search  of  the  parameter  regime 
identified  by  Smith  and  Burggraf  (1985)  no  solutions  could  be  found,  but  this  does  not  mean 
that  they  do  not  exist. 

The  next  calculation  required  concerned  the  constant  Dj  defined  by  (2.20).  In  order  to 
calculate  flj  from  (2.20)  it  is  necessary  to  compute  the  adjoint  function  J  =  (M,  N,  S,  T).  In 
fact  it  is  easier  to  solve  the  problem  for  («i,vi,j»i,  Ai)  directly  and  find  the  value  of  flj  which 
enables  all  the  required  boundary  conditions  to  be  satisfied.  It  was  easier  to  compute  fij  in  this 
way  because  the  system  for  can  be  solved  using  essentially  the  same  iteration 

method  as  used  above  for  the  solution  of  (2.12).  In  Figure  6  we  show  the  dependence  of  K  on 
the  wavenumber  a.  The  fact  that  K  is  singular  at  «  =  ac  =  1.0145  is  a  direct  consequence  of 
the  fact  that  fl'(oc)  =  0.  In  Figure  7  we  show  the  dependence  of  Ao  on  a  and  we  observe  that 
Aoo  is  respectively  negative  and  positive  on  the  lower  and  upper  branches  of  Figure  3.  Here 
the  upper  and  lower  branches  correspond  to  points  on  Figure  3  which  are  respectively  above 
or  below  E.  The  singularity  in  Aoa  is  due  to  the  fact  that  Ao  continues  to  decrease  when  a 
passes  through  oc.  The  fact  that  both  K  and  A0a  change  sign  at  ac  means  that  viscous  effects 
have  essentially  the  same  destabilizing  role  on  the  upper  and  lower  branches  when  uniform 
wavetrains  are  considered;  see  the  following  section. 

Now  let  us  discuss  the  implications  of  our  calculations  for  the  evolution  equation  (2.22) 
which  we  recall  determines  the  wavenumber  a  correct  upto  order  The  term  on  the  right 
hand  side  of  (2.22)  is  due  to  viscous  effects  and  the  results  of  Figures  6,  7  imply  that  viscous 
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effects  are  destabilizing.  The  zeroth  order  approximation  to  (2.22)  yields 


da  da 

W  +  ,J,dx=0 


(3.4) 


where  ug  is  the  group  velocity.  We  see  from  Figure  3  that  the  group  velocity  is  negative  for 
the  upper  branch  and  positive  otherwise.  This  suggests  that  the  upper  branch  solutions  are 
physically  irrelevant  since  their  energy  propagates  upstream.  In  fact,  the  form  of  Figures  2,3, 
and  the  known  result  about  the  stability  at  small  amplitude  of  Tollmien-Schlichting  waves, 
Smith  (1979b),  suggests  that  solutions  corresponding  to  the  upper  and  lower  branches  would 
be  found  to  be  unstable  and  stable  respectively  if  a  Floquet  analysis  of  them  were  carried  out. 

Suppose  then  that  we  consider  the  evolution  of  disturbances  cdrresponding  to  the  lower 
branch  of  Figure  3.  If  at  T  —  0  we  are  given 


a  =  a(X) 


then  for  T  positive  we  have 

a  =  a(X-ug(a)T) 

which  determines  a  implicitly  since  ug  is  a  function  of  a.  It  is  well-known,  eg  Whitham 
(1974),  that  for  positive  u>'g  the  above  solution  will  become  multivalued  after  a  finite  time  if  the 
initial  data  has  a  compressive  part.  This  suggests  that  finite  amplitude  Tollmien-Schlichting 
waves  will  develop  discontinuities  in  wavenumber  and  frequency  as  they  propagate  downstream. 
When  such  shocklike  structures  develop  (3.4)  is  no  longer  valid,  and  the  viscous  term  must  be 
brought  into  play.  We  expect  that  the  situation  then  is  similar  to  that  for  Burgers  equation, 
see  Whitham  (1974),  where  viscous  effects  smooth  out  shocklike  solutions  but  do  not  prevent 
their  development.  However  until  a  numerical  treatment  of  (3.4)  is  carried  out  such  remarks 
should  be  treated  as  speculative.  Now  we  shall  concentrate  on  a  case  where  more  analytical 
progress  is  possible  and  investigate  nearly  uniform  wavetrains. 
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4  Uniform  wavet rains  and  their  stability 


Suppose  that  (ao,flo(<*o))  is  some  point  on  the  curve  shown  in  Figure  3.  The  corresponding 
wave  with 

0  =  ©(qqA'  -  ftoCacOT) 


corresponds  to  a  constant  frequency / wavenumber  solution  of  the  full  2D  triple  deck  problem 
for  Tollmien-Schlichting  waves.  The  stability  of  this  system  can  be  readily  investigated  by  use 
of  the  phase  equation  (2.22).  We  first  write  a  =  ao  -f  A,  where  A  is  small,  and  then  (2.22) 
becomes 

dA  .dA  .A dA  _x,3vt  , dAo  d  fx 

8T  +  o(«o)^  +  D0(«o)A5;Jr  -  K{a0)  ^  &x  J ^  ^  _  xy/3ds 

+0(61'3  A2,  A3). 


Note  that  KAoa  is  negative  on  both  the  lower  and  upper  branches  respectively  of  Figure  3. 

We  can  eliminate  the  term  proportional  ft^ao)  by  an  appropriate  Galilean  transformation. 
If  we  then  take  T  =  0(£~1/,:}),  A  ~  61/3  with  X  =  0(1)  then,  the  limit  6  — ►  0,  a  suitably 
rescaled  version  of  the  above  equation  is 


^ + _  _  jL  r  _S7_ds 
dr+Adt~  dth  (*-01/3 


(4.1) 


Therefore  the  longwave  instability  of  a  uniform  wavetrain  of  Tollmien-Schlichting  waves  is 
governed  by  the  apparently  new  evolution  equation  (4.1). 

Suppose  now  that  at  r  =  0  there  exists  a  small  initial  perturbation  A  =  Ao(£).  The 
linearized  form  of  (4.1)  shows  that  A  evolves  according  to 


A  =  h  Loo  A°^ exp 


ik£  + 


k2r 


r(2/3)(ifc)2/3 


dk , 


where  Aq  is  the  transform  of  the  initial  data  and 

g-irt/3^,-2/3^  fc  >  q, 


{ik) 


-2/3 


e,,/3|A.|-2/3,  £<0 


It  can  be  seen  that  the  ultimate  form  of  the  disturbance  depends  on  Ag(fc).  More  precisely 
we  see  that  when  Ao(£)  is  sufficiently  concentrated  the  solution  will  develop  a  singularity 
at  a  finite  time.  Thus  for  example  an  initial  disturbance  of  Gaussian  form  will  have  AJ$  ~ 
exp(— A:2)  and  a  bounded  sol  ♦ion  will  occur  for  all  time.  However  if  Aj$  ~  exp  |fc|4/3  the 
solution  will  become  unbounded  after  a  finite  time.  This  is  an  important  result  because  it  says 
that  a  constant  wavenumber /frequency  solution  of  the  two-dimensional  triple-deck  problem  for 
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Tollmien-Schlichting  waves  is  always  unstable  to  a  long  wave  instability.  This  is  not  uncommon 
in  physical  problems,  eg  the  Stokes  water  wave,  nonlinear  optics.  For  a  discussion  of  such 
modulational  instabilities  the  reader  is  referred  to  Whitham  (1974). 

A  possible  structure  for  A  when  this  singularity  is  encountered  is  obtained  by  writing 

A  =  (r-r)»*(«,  ^=1^75 

where  the  singularity  occurs  when  r  — ♦  7 _  and  is  localized  around  £  =  £.  The  constant  A  will 
depend  on  the  initial  conditions,  but  we  expect  that  only  solutions  with  A  >  0  w:  possible. 

The  function  ^  is  then  found  to  satisfy 

1T  3 (t><M  d  *'(0)  „ 

+  4  d<f>  ~  d^h  3<l0'  (  '2) 

We  seek  a  continuous  solution  of  (4.2)  with  a  constant,  for  4>  >  0  and  allow  to  grow 

algebraically  when  (f>  — ►  — oo.  It  is  then  convenient  to  write  a  —  —<f>  and  let 


¥(<£)  =  #(-<r)  =  ¥o  {/(*)+!} 


so  that  f(<r)  satisfies 


a  </(„>+. }-?/v)=££ 


f(t)dt 


,  -oo  <  a  <  oo, 


4'  d<r  J-oo(<T-ty/3’ - - 

and  then  f(cr)  =  0,  a  <  0,  whilst  for  a  positive  /  is  determined  by 

„{/(*)- i-»/V))  =  f  («) 

When  n  =  4A/3  =  4/3  this  equation  has  an  exact  solution  proportional  to  a 4/3.  If  this  solution 
is  written  in  terms  of  the  original  variables  we  obtain 

f  *o,  <£>0, 


l  #o  +  (-<t>)4/3,  d><0, 

For  other  values  of  p  we  can  find  a  solution  of  (4.3)  by  taking  a  Laplace  transform.  After  some 
calculation  we  find  that  /(<r)  may  then  be  written  in  the  form 


where 


f(n,  _  ±_  r+,o°  q(p)  , 

27TI  Jy-ioo  / (p) 

I(p )  =  qO-'e-WW’dq. 

J  V 


It  is  then  easy  to  show  from  (4.5)  that  when  a  — *  oo 


f(cT)  ~ 
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Thus  in  terms  of  the  original  variables  we  obtain 


~  (-<£)  3  ,  <(>  — »  -oo. 

which  incidentally  confirms  the  exact  solution  (4.4).  Thus  we  have  constructed  solutions  of 

(4.2)  for  which  —  is  singular  with  0  <  A  <  3/4  whilst  the  solution  for  A  =  0  corresponds  to 
d<p 

a  finite  jump  in  between  (f>  =  ±oo  so  that  the  wavenumber  perturbation  (and  therefore  the 
associated  velocity  displacement  function)  has  a  shock  structure. 

In  the  absence  of  a  full  numerical  investigation  of  the  evolution  equation  (4.1)  we  cannot 
say  which  of  the  singular  solutions  will  be  excited.  Indeed  nothing  in  our  discussion  so  far 
has  ruled  out  the  possibility  of  finite  amplitude  equilibrium  solutions  of  (4.1).  However,  if  we 
multiply  (4.1)  by  X  and  integrate  from  £  =  -oo  to  £  =  -foo,  and  use  Plancherets  formula  to 
simplify  the  contribution  from  the  right  hand  side  of  (4.1),  we  obtain 

V5 


(4.6) 


where  X *  is  the  Fourier  transform  of  X.  It  follows  from  (4.6)  that  equilibrium  solutions 
(constant  or  in  the  form  of  traveling  waves)  are  ruled  out. 

Breakdown  of  the  full  nonlinear  problem 

Now  let  us  consider  a  possible  breakdown  form  for  the  full  nonlinear  system  (4.1).  If  the 
breakdown  is  governed  by  the  inviscid  form  of  the  equation  then,  following  Brotherton-Ratcliffe 
and  Smith  (1987),  we  can,  after  a  suitable  shift  of  origin,  write 


A  =  |r 


n— 1 


Qq(tj)  + 


with  Tj  =  -  ~r  .  The  term  on  the  right  hand  side  of  (4.1)  is  then  negligible  for  n  <  3/4  and  Qq 


is  then  given  implicitly  by 


V  =  ~Qo  -  eoQo> 


(4.7) 


with  eo  a  positive  constant.  However  (4.7)  only  determines  Q0  as  a  single-valued  function  of  7/ 

when  n  =  — — — ,  L  =  3, 5, 7,  •  •  •  so  that  this  type  of  structure  is  not  possible.  However  we  can 
L/  —  1 

take  n  =  3/4  in  which  case  Q o  satisfies 


(4.8) 


The  above  integral  equation  must  then  be  solved  numerically.  We  postpone  that  numerical 
investigation  until  we  have  carried  out  full  numerical  simulations  of  (4.1);  such  an  investigation 
will  be  reported  on  in  due  course. 
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5  A  direct  evaluation  of  the  wavenumber  modulation  equa¬ 
tion 

We  shall  now  use  a  multiple  scale  approach  to  derive  (4.1)  directly  from  the  two-dimensional 
triple  deck  equations  (2.1).  Again  the  major  effect  of  the  modulation  is  to  introduce  a  layer 

t 

of  depth  <5-1/3  sitting  on  top  of  the  lower  deck.  We  again  define  slow  variables  X  and  T  by 
writing 


X  =  6x, 
T  =  St. 


It  is  then  convenient  to  define  X,  Y  by 


X  =  X  -  UgT, 
T  =  6X'3T, 


(5.1a) 

(5.1b) 


(5.2a) 

(5.2b) 


where  ug  is  a  group  velocity  to  be  determined  at  higher  order  in  our  expansion  procedure. 
Suppose  then  that  we  seek  a  solution  of  the  triple-deck  equations  which  is  periodic  in  (f>  = 
ax  —  Sit  where  a  and  SI  are  now  taken  to  be  constant.  In  the  lower  deck  we  expand  u  in  the 
form 

«  =  £  an6n'3Bn(X,t)dnU°{*'y)  +  64'3  [u4(cf>,y,X,f)\ 


n=0 


d<t>n 

+  65'3u5(<t>,y,X,f)  + 


(5.3) 


together  with  similar  expansions  for  v  and  p.  We  note  at  this  stage  that  the  summation  term 
in  (5.3)  arises  because  of  the  translational  invariance  of  a  solution  of  the  two-dimensional 
triple-deck  equations.  In  addition  we  note  that  the  first  correction  to  the  underlying  mean 
state  dependent  on  X  arises  at  0(6 4/f3)  in  (5.3).  It  should  also  be  stressed  at  this  stage  that 
uo,  the  0(1)  term  in  (5.3),  is  independent  of  the  slow  scales  X  and  T  and  that  B(X,T)  is  an 
amplitude  function  to  be  determined.  Thus  in  (5.3)  we  identify  the  term  SX^3B—^-  as  a  small 


du0 


d<f> 


amplitude  perturbation  to  the  periodic  flow  uo (<f>,  y)-  The  eigenfunction  -r—  occurs  because 

o<p 

of  the  translational  invariance  of  any  ^-periodic  solution  of  the  triple  deck  problem.  For  our 
purposes  here  it  is  sufficient  for  us  to  consider  the  partial  differential  equations  to  determine 
(«4,V4,p4)  and  (u5,V5,p5).  If  the  expansions  for  u,v,p  are  substituted  into  the  triple  deck 
equations,  and  the  appropriate  change  of  variables  made,  then  we  find  that  («4,  u4,  P4)  satisfies 


du4  dv4 


(5.4a) 


d2u4  dp4  du4  du0 

~dtf  ~  a~dt  ~  au°!ty  ~  aU4~dt 


V4lfy  ~  V°lfy  =  ^-u;sli0  +  au°u°<i]  Bx  ■  (5-4b) 
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These  equations  must  be  solved  subject  to  U4  =  =  0,  y  =  0  whilst  for  large  y  the  appropriate 

conditions  are 


u4  ->  A(<f>,X,T), 


(5.5a) 


with 


1  r 

P4  =  -  f 

air  J-c 


dA 

■57 


( <t>~8 ) 


ds 


(5.6a) 


Since  the  homogeneous  form  of  the  system  for  has  the  solution  (U4 ,v4,p4)  = 

d 

— («o,  vo,po)  it  follows  that  a  solution  exists  only  if  an  orthogonality  condition  is  satisfied. 
od> 

This  may  be  written  down  following  the  procedure  used  in  §4,  it  is  sufficient  here  to  note  that 
the  condition  determines  the  group  velocity  u>g  and  that  the  expansion  obtained  is  identical  to 
that  which  it  was  derived  in  §3  after  making  a  perturbation  in  the  wavenumber.  The  solution 
of  (5.4)  is  clearly  of  the  form 


(U4,t?4,P4)  —  (1/4,  V4,  ^*4),  ^4  —  BfA, 

where  (U4,  V4,  P4)  and  A  are  independent  of  X,  T.  The  other  main  feature  to  appreciate  about 
the  solution  of  the  0(64/3)  problem  is  that  at  the  edge  of  the  y  —  0(1)  region  U4  may  be 
expressed  in  the  form 

U4~  A  =  (Am  +  Ap(4>))Bj( 

where  Am  is  independent  of  <j)  and  therefore  corresponds  to  a  mean  flow  correction.  The 
reduction  to  zero  of  Am  is  achieved  in  the  outer  0(S~1^3)  region  in  a  similar  manner  to  that 
found  in  §3.  We  define  the  variable  tj  by 

t)  =  S1,3y 


and  in  the  outer  6  a/3  layer  u  is  expanded  in  the  form 


«  =  7^75  +  £  an6n'3Bn(X,f  &(4>,  00)  +  64'3UM(X,v )  +  *5/3^  + 


tfl/3 


»=1 


d<t>- 


where  Um  is  the  mean  flow  correction  driven  by  Am-  The  mean  flow  correction  in  the  ij 
direction  is  62Vm  and  the  linear  problem  to  determine  ( Um »  Vm)  is  found  to  be 


dUM 
OX 
d2UM 
di i2 
Um 


dVM 


-  V- 


dr/ 

8Um 


=  0, 


-Vm  =  0, 


ex 

=  AmB%,Vm  =  0,tj=  0 


Um  -*•  0,t?->oo. 


(5.7) 


The  system  (5.7)  is  identical  to  that  obtained  in  Section  2  and  therefore  may  be  solved  again 
using  a  Fourier  Transform  technique,  the  solution  is  not  repeated  here.  The  mean  flow  at  order 
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6 4/3  then  interacts  with  the  0(£°)  flow  to  produce  an  0(6S/3)  correction  to  the  outer  boundary 
condition  for  the  disturbed  flow  in  the  y  =  0(1)  layer.  Again  the  analysis  follows  closely  that 
of  §2  so  we  do  not  repeat  it  here,  we  find  that  {us,v5,p5)  must  satisfy 


dus 

d<f> 


+ 


r 

dX  Jx 


%(«.*) 

( s  —  i)1/3 


ds, 


where  A  is  a  constant.  In  the  y  =  0(1)  region  (1*5,175,  ps)  is  found  to  satisfy  (5.4a, b)  but  with 
the  right  hand  sides  of  these  equations  replaced  by  [•],  and 

„  ,  „  ,  r.  n  f  dUQ  du4  duo  du4  ) 

[•]  +  Bfuo<,  +  BBk  |t;4_  +  Uo_  +  V4—  +  u0—  +  •••} 

respectively.  Here  [•]  denotes  terms  which  don’t  contribute  to  the  solvability  condition  which 
the  system  for  (u5,  v5,p5)  must  satisfy.  The  required  condition  is 


8B  dB  t  3  /“  (S,f)  , 

dr  ax  ax  Jx  (s-xy/ 3 


(5.8) 


where  ,  A  are  constants  and  a  suitable  change  of  variables  enables  us  to  recover  (4.1).  Thus 
T  in  (4.1)  can  be  interpreted  as  being  either  a  wavenumber  or  amplitude  perturbation  and  the 
two  approaches  are  shown  to  be  consistent. 


6  The  phase  equation  approach  at  finite  Reynolds  numbers 


In  §2  we  derived  the  phase  equation  for  Tollmien-Schlichting  waves  in  a  Blasius  boundary  layer. 

Such  a  boundary  layer  exists  only  at  asymptotically  large  values  of  the  Reynolds  number  and 

it  was  therefore  appropriate  to  utilize  the  largeness  of  the  Reynolds  numbers  in  the  description 

of  the  instability  wave.  However  linear  descriptions  of  the  evolution  of  Tollmien-Schlichting 

waves  at  finite  Reynolds  numbers  have  been  given  by  for  example  Gaster  (1974).  Though  such 

approaches  do  not  give  formal  asymptotic  approximations  to  the  equations  of  motion  it  appears 

that  they  correctly  predict  the  essential  physics  of  the  linear  growth  of  Tollmien-Schlichting 

waves.  Similarly  large  scale  numerical  simulations  of  nonlinear  growth  of  Tollmien-Schlichting 

waves  at  finite  Reynolds  numbers  have  proved  equally  successful  at  reproducing  experimental 

results,  e.g  Wray  and  Hussaini  (1984).  Here  we  wish  to  investigate  the  phase  equation  approach 

at  finite  Reynolds  numbers  but,  in  order  to  keep  our  asymptotic  analysis  formally  correct,  we 

choose  to  work  with  a  parallel  boundary  layer  which  is  an  exact  solution  of  the  Navier-Stokes 

equations  at  all  Reynolds  numbers.  We  refer  to  the  asymptotic  suction  boundary  layer  which 

has  been  investigated  in  the  weakly  nonlinear  regime  by  Hocking  (1975).  Suppose  then  that 

the  freestream  speed  is  (Jo  and  the  suction  velocity  is  Vo-  We  define  a  reference  length  L  =  — - 

ko 

and  define  the  Reynolds  number 

R=Uo 
R  Vo 
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but  we  assume  that  R  =  0(1)  in  this  section.  Since  we  restrict  ourselves  to  two-dimensional 
disturbances  it  is  convenient  for  us  to  define  a  stream  function  0  and  work  with  the  vorticity 
equation  in  the  form 

^  +  ejp^  =  L 

at  d(x,y)  R 


which  is  to  be  solved  subject  to 


i’y  =  o,  =  i>  y  =  o, 


(6.2a) 


i’y  -»  1,  y  —  oc 

In  the  absence  of  an  instability  the  stream  function  V’  is  given  by 


(6.2b) 


V>  =  Mv)  =  y  +  e  y  -  x.  (6.3) 

This  flow  is  unstable  to  two-dimensional  Tollmien-Schlichting  waves  for  R  >  54370.0  and  the 
band  of  unstable  wavenumbers  tends  to  zero  when  R  — ►  oo.  Here  we  assume  R  is  0(1)  and 
assume  that  an  0(1)  amplitude  wavesystem  is  superimposed  on  (6.3).  At  finite  R  the  mean  flow 
driven  by  the  wave  system  is  confined  to  the  boundary  layer.  Therefore  no  outer  adjustment 
layer  is  required  even  when  the  wavesystem  evolves  slowly  in  the  downstream  direction,  see 
Hocking  (1975)  for  a  discussion  of  this  point.  Suppose  that  X,T  are  defined  by 


then  we  define 


and  expand  il  in  the  form 


X  —  Sx,  T  =  St, 


Q  =  -0(X,T),  a  =  0x,  il=  -0T, 


fl  —  fig  4"  4" 


We  seek  a  solution  of  (6.1)  by  writing 


i>  =  MX,  y ,  ©,  T)  +  SMX,  y,  0,  T)  4- 


and  the  leading  problem  is  determined  by  solving 


r,  1^4, 

~n°-ae~  +  a  mV)  =R  '*°' 

M  =  0,  «V>oe  =  1,  y  =  0, 

M  -» l,  y  ->  oo, 


_2  _  d 2  2  d 2 

Vl  =  %2  +Q  ae2’ 
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and  periodicity  in  0.  The  required  frequency  and  wavenumber,  fio  and  a,  must  be  found 
numerically  once  some  measure  of  the  disturbance  size  is  specified.  We  will  not  attempt  such 
a  calculation  here  but  we  note  from  the  finite  Reynolds  number  calculation  of  Hocking  (1975) 
and  the  large  Reynolds  number  theory  of  Smith  (1979a,b)  that  both  sub  and  supercritical 
bifurcations  to  finite  amplitude  Tollmien-Schlichting  waves  are  possible.  For  the  purposes  of 
our  discussion  we  simply  assume  that  the  nonlinear  eigenrelation  fio  =  flo(<*,  R)  is  known.  At 
order  6  we  find  that  V'i  is  determined  by 


O  dViV>i  _  d(xi>o,v2M 

°  dQ  °  d(G,y )  a  d(G,y) 


<6-6' 


where 


M  = 


_  d(V>0,  Vfyo)  2  2 date)  , 

d(a,y)  d(G,y)  R  dadG 


V2 


^o- 


The  system  must  be  solved  subject  to  periodicity  in  0  whilst  the  boundary  conditions  in  y 


are 


dib\ 

V>i  =  =0,  y  =  0 

i  —  q(X),  y  —*  oo. 


(6.7a) 

(6.7b) 


Here  q(X)  represents  a  mean  flow  normal  to  the  wall  at  infinity.  This  flow  is  essentially 
driven  through  the  equation  of  continuity  by  the  0(1)  streamwise  velocity  component.  A  solv¬ 
ability  condition  is  required  if  (6.6)-(6.7)  is  to  have  a  solution  since  the  translational  invariance 

of  (6.5)  means  that  V>i  =  is  a  solution  of  the  homogeneous  form  of  (6.6)-(6.7).  It  is  worth 

a0 

pointing  out  at  this  stage  that,  if  we  were  performing  a  calculation  in  a  region  of  finite  depth, 
then  a  pressure  eigenfunction  would  have  to  be  allowed  for  at  leading  order  in  order  that  q 
should  be  reduced  to  zero.  The  requirement  for  such  an  eigenfunction  is  well-known  in  weakly 
nonlinear  stability  theory;  see  for  example  Davey,  Hocking  and  Stewartson  (1974)  or  DiPrima 
and  Stuart  (1975).  The  solvability  condition  can  be  found  by  writing 


in  which  case  the  homogeneous  form  of  (6.6)  is 

°^2+|;fi2+c'2=0' 

where  A,  B  and  C  are  6x6  matrices  defined  by 

'  1  0  0  0  0  0  ' 

0  0  0  0  0  0 

_  0  1  0  0  0  0 

0  0  0  1  0  0  ’ 

0  0  0  0  0  0 

0  0  0  0  1  0 
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0  0  0  0  0  0 

1  0  0  0  0  0 

_  0  0  1  0  0  0 

0  0  0  0  0  0 

0  0  0  1  0  0 

.0  0  0  0  0  1 

0-1  000 

00  -100 

0  0  0  -1  0 

0  0  0  0  -1 

0  0  0  0  0 

0  oy  a#V?0oe  0  R[aip0y  -  SI] 

The  system  adjoint  to  that  given  above  is 

-°lkAT<}-TyBTCi  +  CTq  =  0’ 

together  with  conditions  of  periodicity  in  0  and 

<75  =  9e,  V  =  0,  oo. 

The  condition  that  (6.6)-(6.7)  has  a  solution  then  becomes 

n>  =  -m'" 

with 

C 

The  phase  condition  Six  +  =  0  correct  upto  0(6 )  becomes 

«T  +  ax fi0(«)  =  Tpj£  {ox  •  h  («)}  6. 

In  order  to  examine  the  stability  of  a  uniform  wavestream  solution  we  write 

a  =  a0  +  A 

with  A  <<  «o,  and  (6.11)  after  a  suitable  change  of  scales  becomes 

At  +  AA^  =  ±A^. 


0 

0 

0 

0 

-1 

-Raipo& 


(6.8) 

(6.9) 

(6.10) 

(6.11) 

(6.12) 


Thus  we  obtain  the  surprising  result  that  the  modulational  instability  of  a  two-dimensional 
wavesystem  in  a  boundary  layer  at  finite  Reynolds  numbers  is  governed  by  Burgers  equation. 
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In  (6.12)  the  ±  signs  correspond  to  the  cases  when  diffusion  effects  are  stabilizing/destabilizing 
respectively.  Without  calculating  fti  we  do  not  know  which  sign  is  appropriate  for  the  problem 
under  consideration  here  so  we  will  discuss  both  possibilities.  However  we  can  simply  quote 
the  known  results  about  Burgers  equation  for  each  case  and  the  implications  for  the  stability 
of  a  uniform  wavetrain  are  essentially  the  same.  A  full  discussion  of  results  quoted  below  can 
be  found  in  Whitham  (1974)  and  Howard  and  Kopell  (1977). 

If  the  positive  sign  is  taken  the  solution  remains  bounded  for  all  time  and  indeed  localized 
or  periodic  solution  of  (6.12)  tend  to  zero  when  T  increases.  However  even  in  the  diffusively 
stable  case  weak  shock  solutions  of  (6.12)  can  develop;  a  discussion  of  this  possibility  is  given 
by  Whitham  (1974),  Howard  and  Kopell  (1977).  Thus  for  any  gi^en  uniform  wavetrain  whose 
instability  is  governed  by  (6.12)  with  the  positive  sign,  an  initial  disturbance  can  be  found 
which  does  not  decay  to  zero  at  large  times.  The  uniform  wavetrain  is  therefore  modulationally 
unstable. 

If  the  negative  sign  is  taken  in  (6.12)  viscous  effects  are  destabilizing  and  in  fact  finite  time 
singularities  are  developed  from  a  broad  range  of  initial  conditions.  Again  it  follows  that  the 
uniform  wavetrain  is  unstable. 

We  conclude  that  at  finite  Reynolds  numbers  modulational  effects  will  either  cause  a  finite 
time  singularity  to  develop  if  viscous  effects  are  destabilizing,  or  cause  shock  discontinuities 
in  wavenumber  and  frequency  in  the  stable  case.  In  either  case  uniform  wavetrains  of  two- 
dimensional  Tollmien-Sclichting  waves  at  finite  amplitude  should  not  be  observed  according  to 
our  theory. 

7  Conclusion 

We  have  used  a  phase  equation  approach  to  determine  the  evolution  of  Tollmien-Schlichting 
waves  at  large  and  finite  Reynolds  numbers.  In  the  high  Reynolds  number  case  we  found  that 
finite  amplitude  disturbances,  periodic  in  time  and  space,  apparently  exist  only  for  a  small 
range  of  values  of  the  wavenumber  a.  The  upper  branch  in  Figure  3  describes  modes  with 
negative  group  velocity  and  so  they  are  therefore  of  no  physical  interest.  The  lower  branch  on 
the  other  hand  corresponds  to  waves  with  positive  group  velocity  and  for  these  modes  the  rate 
of  change  of  the  group  velocity  with  a  is  positive.  This  means  that  (2.22),  the  leading  order 
approximation  to  the  phase  equation,  will  develop  discontinuities  after  a  finite  time  for  many 
initial  disturbances.  It  might  be  anticipated  that  viscous  effects,  which  appear  on  the  right 
hand  side  of  (2.22),  might  smooth  out  such  discontinuities.  We  cannot  be  sure  if  that  is  the 
case  until  a  numerical  investigation  of  (2.22)  is  carried  out.  However,  in  §4  we  investigated  the 
particular  case  of  uniform  wavetrains  and  found  that  there  viscous  effects  were  destabilizing 
and  it  seems  likely  that  this  is  also  the  case  for  (2.22).  In  §4  we  saw  that  the  linearized  form 
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of  (4.1)  leads  to  a  finite  time  singularity  of  the  type  found  by  Brotherton-Ratcliffe  and  Smith 
(1987).  For  the  full  nonlinear  problem  (4.6)  shows  that  equilibrium  solutions  are  not  possible 
and  that  _  n  initial  perturbation  cannot  decay  to  zero  at  large  times.  In  effect  this  means  that 
periodic  solutions  of  the  triple-deck  equations  are  modulationally  unstable.  A  possible  form  for 
the  structure  of  a  singular  solution  of  the  breakdown  of  the  nonlinear  form  of  (4.1)  was  found. 
The  structure  found  was  again  based  on  the  structure  found  by  Brotherton-Ratcliffe  and  Smith 
(1987).  The  validity  of  the  structure  must  remain  open  until  a  full  numerical  investigation  of 
(4.1)  has  been  carried  out.  The  question  of  how  the  Navier-Stokes  equations  alter  their  large 
Reynolds  number  structure  in  order  to  remove  the  singularities  of  (4.1)  also  remains  open. 

At  finite  Reynolds  numbers  we  found  that  the  evolution  equation  for  a  periodic  wavetrain 
satisfies  Burgers  equation.  Without  extensive  calculations  we  cannot  say  whether  the  viscous 
term  in  (6. 12)  has  positive  a  or  negative  sign.  If  it  turns  out  to  be  negative  then  viscous  effects 
are  again  destabilizing  and  finite  time  singularities  will  occur.  If  the  sign  is  positive  then  viscous 
effects  are  stabilizing  .  However  (6.12)  can  be  solved  exactly  by  the  Cole-Hopf  transformation 
and  it  is  known,  Whitham  (1974),  that  even  in  the  stable  case  shocks  will  in  general  develop. 
Thus  we  conclude  that  at  large  or  finite  Reynolds  numbers  a  uniform  wavtrain  of  Tollmien- 
Schlichting  waves  will  break  down  with  a  singularity  or  shock  developing  after  a  finite  time. 
This  casts  some  doubt  on  the  validity  of  large-scale  simulations  of  Tollmien-Schlichting  waves 
using  Fourier  series  expansions  in  the  streamwise  direction. 

In  view  of  the  fact  that  our  analysis  has  been  restricted  to  the  two-dimensional  case  it 
is  possible  that  three-dimensional  effects  might  prevent  the  above  predictions  from  occurring 
in  practice.  Nevertheless,  experimental  observations  where  the  Tollmien-Schlichting  wave  is 
driven  by  a  wavemaker  suggest  that  the  first  step  in  the  transition  process  is  the  linear  growth 
of  two-dimensional  Tollmien-Schlichting  waves  followed  by  nonlinear  saturation  and  three- 
dimensional  effects  coming  into  play,  see  for  example  Klebanoff  et  al  (1962). 

As  a  first  step  towards  looking  at  the  effect  of  three-dimensionality  let  us  determine  how 
the  evolution  equations  (4.1),  (6.12)  are  modified  if  a  weak  spanwise  dependence  of  the  wave 
is  allowed  for.  We  shall  determine  the  modifications  using  the  multiple-scale  approach  of  §5 
rather  than  modify  the  phase  equation  formulation.  Since  §5  concerned  Tollmien-Schlichting 
waves  at  large  Reynolds  number  let  us  give  the  essential  details  of  how  the  corresponding 
expansion  procedure  is  formulated  at  finite  Reynolds  numbers.  In  the  absence  of  spanwise 
modulation  the  streamwise  velocity  component  is  written  in  the  form 


u  =  £antfnA  n(X,T) 
n=0 


-  dnu0{<t>,y) 


d<j>n 


+  62u-2(<t>,  y,  X,  T )  +  63u3(<j>,  y,X,T)  + 


(7.1) 


Here  X  =  6(x  —  u>gt)  is  a  slow  variable  obtained  by  moving  to  a  coordinate  system  moving 
downstream  with  the  group  velocity  and  <f>  =  ax  —  Slot.  The  summation  term  arises  because  of 
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the  translational  invariance  of  a  periodic  solution  in  the  streamwise  direction.  At  order  62  the 
group  velocity  is  determined  as  a  solvability  condition  and  u2  will  be  proportional  to  A^.  At 
order  63  terms  proportional  to  AA^  and  a  slow  time  derivative  of  A  will  be  generated.  Viscous 
effects  also  come  into  play  to  generate  terms  proportional  to  A^  and  the  required  solvability 
condition  leads  us  to  (6.12).  In  the  presence  of  spanwise  modulations  the  procedure  described 
above  is  modified  in  the  following  way.  Firstly  we  define  a  slow  variable 


£  =  6z 

where  z  is  the  spanwise  variable  and  allow  the  amplitude  function  A  to  depend  on  £.  The  above 

* 

expansion  for  u  remains  unchanged  but  the  0(6)  term  in  the  corresponding  pressure  expansion 
drives  a  spanwise  flow  of  order  6 2  proportional  to  A^.  This  spanwise  velocity  component  then 
generates  a  term  proportional  to  A^  in  the  0(63)  continuity  equation.  Similar  terms  are 
generated  by  spanwise  diffusion  in  the  x,  y  momentum  equations.  The  solvability  condition 
then  is  modified  to  give 

At  +  AA^  ±A^  +  AiA^.  (7-2) 

Here  Ai  is  a  constant  whose  sign  cannot  be  calculated  without  solving  the  0(1)  problem 
numerically.  A  similar  analysis  for  the  triple-deck  case  gives 


d\  Kd\ 

87  +  Aa« 


oa  _  _d_  _ I 

^  "  df  Jt  (s  -  o1/ 3 


ds  +  A2A££. 


with  A2  a  constant  as  the  required  generalization  of  (4.1).  As  a  starting  point  to  study  the  role 
of  three-dimensional  effects  in  the  two-dimensional  breakdown  structure  w~  plan  to  carry  out 
numerical  simulations  of  (7.2),  (7.3). 

Finally  we  close  with  a  brief  outline  of  how  (2.22),  (4.1)  are  modified  by  nonparallel  flow 
effects.  Within  the  framework  of  (2.1)  nonparallel  effects  manifest  themselves  through  the 
pressure  displacement  law  which  allows  for  elliptic  effects  in  the  streamwise  direction.  However 
the  slow  spatial  evolution  of  the  unperturbed  shear  flow  does  not  enter  (2.1)  since  it  can  be 
scaled  out  of  the  problem.  For  the  weakly  nonlinear  growth  of  Tollmien-Schlichting  waves 
Hall  and  Smith  (1984)  showed  that  nonparallel  effects  are  important  when  the  disturbance 
amplitude  is  0(R~7^32)  and  lead  to  an  amplitude  equation  of  the  form 

^  =  XC-(l+i6)C\Cf,  (7.4) 

with  6  a  real  constant.  Thus  nonparallel  effects  lead  to  the  term  XC  in  (7.4)  thereby  causing 
the  increased  linear  exponential  growth  of  a  small  disturbance  as  it  evolves  in  X .  Any  small 
disturbance  amplifying  by  this  means  eventually  becomes  nonlinear  and  for  large  X  has  |(7|2  ~ 
X .  Let  us  now  show  how  related  terms  can  be  incorporated  into  (2.22)  and  (4.1). 
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We  note  that  (2.11a)  proceeds  in  powers  of  Si  and  that  the  Reynolds  number  has  been 
effectively  scaled  out  of  the  problem  by  our  assumption  that  the  Tollmien-Schlichting  wavesys- 
tem  is  described  by  the  triple-deck  system  (2.1).This  assumption  means  that  the  analysis  given 
so  far  in  this  paper  is  formally  valid  for  £  large  compared  to  any  positive  power  of  c  =  R~i .  In 
order  to  reveal  the  effect  of  boundary  layer  growth  we  now  relax  this  condition  and  see  which 
new  terms  now  play  a  role  in  (2.11a).  Clearly  (2.11a)  must  include  terms  proportional  to 
powers  of  e  because  of  the  higher  order  terms  in  the  triple  deck  expansion  but  in  addition  there 
will  be  a  term  proportional  to  (3S~iX  obtained  by  expanding  the  streaniwise  dependence  of 
the  unperturbed  flow  in  a  Taylor  series  in  the  streamwise  direction.  We  therefore  now  expand 
the  frequency  in  the  form 

ft  =  ft0  +  ^1/3fti  +  eft-2  +  e'2ft3  +  (3S~iX(l4  +  0((3,(66~2 ,6i ).  (7.5) 

The  ordering  of  the  terms  in  the  above  expansion  depends  on  the  relationship  between  S  and  (. 

The  term  proportional  to  ft4  is  the  first  one  dependent  on  the  nonparallel  nature  of  the  basic 

flow.  The  first  significant  distinguished  limit  arises  when  6  ~  c  when  the  terms  proportional  to 

2 

ft3,ft4  become  comparable  but  still  small  compared  to  the  0{6 3)  term.  The  next  significant 

9  2 

stage  is  when  6  decreases  to  c*  in  which  case  the  nonparallel  term  and  the  0(6i )  term  are 

9 

comparable.  However  the  crucial  stage  arises  when  6  decreases  to  e*  in  which  case  the  terms 
proportional  to  fti ,  ft4  play  an  equal  role  and  viscous  and  nonparallel  effects  are  comparable. 
Thus  if  we  write 

e3  =  h^Si , 

with  h\  an  0(1)  constant  then  (2.2),  correct  to  second  order,  now  gives 

da  dSl0da  ,1/3,0fti  dSl4  .  . 

Wf  +  J^dx--*  {dx+','si,  +  h'x~dx  ")  (76) 

A  further  decrease  in  the  size  of  6  means  that  the  term  proportional  to  fti  should  be  dropped 
and  nonparallel  effects  dominate  the  right  hand  side  of  the  equation  for  a.  This  completes  our 
description  of  how  the  phase  equation  changes  its  structure  when  6  is  decreased.  Note  here 
that  when  we  decrease  6  we  are  in  effect  moving  further  away  from  the  initial  location  x *  =  x * 
so  the  above  different  orderings  of  the  right  hand  sides  of  (7.5)  correspond  to  moving  further 
downstream.  Thus  the  final  form  has  the  right  hand  side  of  the  phase  equation  dominated  by 
nonparallel  effects. 

A  similar  procedure  can  be  used  to  determine  the  appropriate  modifications  to  (4.1)  which 
is  the  evolution  equation  for  the  perturbed  wavenumber  of  an  initially  uniform  wavetrain.  Here 
the  crucial  scaling  brings  in  nonparallel  effects  at  the  same  stage  as  the  integral  term  driven 
by  viscous  terms.  The  appropriate  scaling  now  has 

=  ho 
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with  h2  an  0(1)  constant.  Equation  ^4.1)  then  becomes 


dA  dA 

—  +  (A+h3T)  — 


d_  r°° 

dUi 


dA 

ds 

(s-zy/3 


ds  -t-  /14. 


(7.7) 


Here  /i3,/i4  are  constants  proportional  to  hi,  and  /14  can  be  set  equal  to  zero  by  a  change  of 
dependent  variable.  Therefore  the  nonparallel  modulational  equation  for  Tollmien-Schlichting 
waves  in  a  growing  boundary  layer  is  (7.6)  with  /14  =  0,  the  numerical  solution  of  that  equation 
will  be  reported  on  in  due  course. 
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FIGURES 


Figure  1  The  dependence  of  a  on  \Jcq  +  <%. 

Figure  2  The  dependence  of  flo  on  \Jcl+ ijj. 

Figure  3  The  dependence  of  fio  on  a.  The  symbols  denote  the  results  of  Conlisk  et  al. 

Figure  4  The  shear  stress  as  a  function  of  0  for  flo  =  2.2995,2.3041,2.3125,2.3245,2.3398, 
2.3575, 2.3763, 2.395, 2.4124, 2.4275 

Figure  5  The  shear  stress  44^(0)  as  a  function  of  flo  for  the  first  sixteeen  modes. 

Figure  6  The  dependence  of  A”  on  a. 

Figure  7  The  dependence  of  Aq  on  a. 
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